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Shor's factorisation algorithm is a combination of classical pre- and post-processing and a quantum 
period finding (QPF) subroutine which allows an exponential speed up over classical factoring algo- 
rithms. We consider the stability of this subroutine when exposed to a discrete error model that acts 
to perturb the computational trajectory of a quantum computer. Through detailed state vector sim- 
ulations of an appropriate quantum circuit, we show that the error locations within the circuit itself 
heavily influences the probability of success of the QPF subroutine. The results also indicate that the 
naive estimate of required component precision is too conservative. 
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1 Introduction 

The investigation and implementation of large scale quantum algorithms is arguably of enormous im- 
portance to the field of quantum information processing. The seminal work by Shor in 1994 was 
the first example of a complex and large scale algorithm that was able to efficiently solve a classically 
intractable problem. Since Shor's discovery, the construction of a large scale quantum computer (QC) 
has been an area of intense research. Currently there are many different proposals for constructing 
such a device [2 3 j, but despite significant progress, the issue of decoherence and imperfect gate de- 
sign begs the question of whether such a large and complex algorithm can be experimentally realized 
beyond trivial problem sizes. 

The development of quantum error correction (QEC) Q] [5] |6) and fault-tolerant quantum compu- 
tation Q|SJ|9l has shown theoretically how large scale algorithms can be implemented on imperfect 
devices. However, without a working QC, detailed classical simulations of QEC and quantum algo- 
rithms constitute the only method for reliable information regarding the behaviour of such schemes 
and the ease in which they can be implemented on physical systems. The issue of appropriate use of 
QEC and the construction of arbitrary fault-tolerant gates 1 10 1 still requires detailed knowledge of the 
behaviour of the underlying algorithm in order to tailor these schemes appropriately. 

For large scale quantum algorithms, the general method of analysis is to assume that all components 
within the algorithm have a precision of as 1/ n p , where n p = KQ represents the number of locations 



1 



2 Title Robustness ofShor's algorithm. 



where an error can occur during an algorithm utilising Q qubits and K elementary steps (depth of the 
circuit). This estimate implies that a single error anywhere during calculation will result in failure. For 
small quantum circuits, this approximation is not an obstacle in component design. However, for more 
complex circuits, where qubits may be coupled in highly non-trivial ways, it is not obvious that such a 
naive estimate is sufficient. In fact, our results show that they are not. In our analysis we examine the 
quantum period finding (QPF) subroutine, which lies at the heart of Shor's algorithm, in the presence 
of discrete errors. The choice of the QPF subroutine in this analysis is due to its importance to the 
field of quantum computing and because it is a good example of a well known, non-trivial algorithm. 

Quantum circuits to factor large integers, for example a 128-bit number, require n p of the order 
10 7 — 10 10 depending on the specific circuit used. Engineering quantum gates with failure rates 
of 10~ 7 — 10~ 10 is currently far from being experimentally realized in any of the numerous archi- 
tectures currently proposed. Our simulations show that the 1 jn p precision requirement is not strictly 
required. We find evidence for a required precision of P(L)/n p , where P(L) is a monotonically in- 
creasing function of L, the binary length of the composite number, which is at least linear. This slower 
scaling increases the error rate at which quantum processing (as opposed to classical randomness) can 
be observed. 

Several authors have previously examined the effects of errors on Shor's algorithm 11311141 IT51 . These 
simulations are often limited to specific sections of the entire circuit, or to other sources of error such 
as phase drifts on idle qubits, imperfect gate operations or aspects relating to quantum chaos. Chuang 
et al [ 16 1 was one of the first to look at the error stability of Shor's algorithm, analytically, under the 
effects of environmental coupling. Miquel et al | 17 1 examined the stability of Shor's algorithm using 
an identical error model to that used in this investigation. However, the stability of the algorithm was 
only investigated for a single problem size and did not investigate how the stability changes as the 
problem size increases. 

Several architectures, most notably solid state models, are restricted to a single line of qubits with 
nearest neighbour interactions only. The issue of whether the QPF subroutine can be implemented 
on such linear nearest neighbour (LNN) architectures is also investigated and compared with circuits 
designed for architectures that can interact arbitrary pairs of qubits (non-LNN). We find that if LNN 
circuits can be designed with comparable values of n p , the stability will be similar. 

In this paper we examine specific circuits for both LNN and non-LNN architectures in the presence of 
a discrete error model, in order to determine: 

• The degree to which the final required state of the computer is affected by small changes in the 
computational trajectory caused by these errors. 

• The impact of a LNN architecture on the reliability of the QPF subroutine. 

• If the 1 jn p bound for component precision remains absolute for various problem sizes. 

The paper is organised as follows. Section|2]examines the underlying theory behind Shor's algorithm, 
the QPF subroutine and how success is defined. Section|3]details the error model and issues relating to 
simulations. Section|4]present simulation results, examining the stability of the QPF subroutine near 
the l/n p lower bound for both LNN and non-LNN circuits. Finally we present a brief analysis that 
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examines the consequence of various additional scalings of component precision when attempting to 
observe quantum processing for small instances of the QPF subroutine. 

2 Shor's algorithm 

As several papers detail the major steps of Shor's algorithm | fl ll8M19l . we provide an overview for 
the sake of completeness and to introduce notation. We first consider a given composite number 
TV = N1N2 which has a binary length L = log 2 (AT). To factorise this number, we consider the func- 
tion f(k) = x k modN, where k 6 Z and x is a randomly chosen integer such that 1 < x < N and 
gcd(iV, x) = 1 (gcd = greatest common divisor). The QPF subroutine of Shor's algorithm determines 
the period of f(k). i.e. to find the integer r > such that f(r) — 1. This QPF subroutine is the 
quantum component of Shor's algorithm. The complete algorithm is composed of both the QPF sub- 
routine and several pre and post processing operations that can be performed in polynomial time using 
classical techniques. These classical steps, detailed by several authors GJHQjO, can be implemented 
in polynomial time and for our purposes we assume that these steps can be implemented with no error. 
Once the period of f(k) is found, the factors of N can be calculated as Ni — gcd(/(r/2) — 1, N) and 
N2 = gcd(/(r/2) + 1, N), conditional on r being even and f(r/2) ^ N — 1, 

In general, to factorise a number of binary length L, 3L qubits are initialised to the state |0) 2_l |0) z,. 
For clarity we have broken these 3L qubits into 2L qubits to store the values k and L qubits to store the 
function evaluations, x k modN. After initialisation, a Hadamard transform is performed on each of the 
2L qubits, placing the k register into an equal superposition of all binary numbers from — > 2 2L — 1, 

2 2L — 1 

|0) 2L |0) L £ \k) 2L \0)L. (1) 

k=0 

Step three is to apply the function f(k) on the L qubit register, conditional on the values k. The state 
of the computer is transformed to, 

2 2t_ 1 2 2i -l 

¥ £ \k) 2L \0)L \k) 2 L\x k modN) L . (2) 

fe=0 k=0 

The next step is to measure the L qubit register. This step can actually be omitted when implementing 
the algorithm, however we introduce it to show how the period, r, appears within the procedure. After 
measurement the qubit register collapses to, 

2 2L -1 2 2L /r-l 

^ £ \k) 2L \f(k)) L ^^ £ \k +nr) 2L \f ) L . (3) 

k=0 n=0 

Where r is the period of /, f is the measured value and k is the smallest value of k such that 
/o = f{ko)- We now apply a quantum Fourier transform (QFT) to the k register. The state of the 
computer after the application of the QFT becomes, 

( 2 ^_i) (2 2t /r-l) . . 

E E exp ^(^o + nr) b) 2L |/ ) L . (4) 

j=0 n=0 ^ ' 
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If we now measure the k register, we will return a value of j with probability, 

P(j,r,L) 



^2L 2^ eX P(^ nr 



(5) 



Eq. |5]is strongly peaked at certain values of j. If the period r perfectly divides 2 2L then Eq. |5]can 
be evaluated exactly, with the probability of observing j = c2 2L /r for < c < r being 1/r, and if 
j 7^ c2 2L /r [Fig. Ha)]- If r * s not a perfect divisor of 2 2L , then the peaks of Eq. |5]become slightly 
broader, [Fig. CQb)], and classical methods can be utilised in order to determine r from the measured 
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Fig. 1. Plot of Eq.gjfor the case, 2 2L = 256 with a) r = 8 and b) i 
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value of j . Given several measured integer values around these non-integer peaks a continued fractions 
method can be employed to determine r (2]|19I . The probability of success s for Shor's algorithm is 
generally defined as, 



s(L,r)= p(j,L,r). 



(6) 



{useful j] 



{useful j} is the set, j = [c2 2L /rj, j = \c2 2L /r], < c < r, where |_ | ] denote rounding down 
and up respectively and p(j, L, r) is defined via Eq. [5] Using this definition of s we determine the 
period after 0(1/ s) calls to the subroutine. 

Many circuits have been proposed in order to implement the QPF subroutine on a physical quan- 
tum computer, as summarised in tabled Some are optimised for conceptual simplicity |20|, some for 
speed 1211 and some for utilising a minimum number of qubits I 22l l24l . 



This investigation will focus on circuits that require a minimal number of qubits for two reasons. 
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Table 1 . Number of qubits required (Q) and circuit depth (K) of different implementations of the QPF subroutine. 
Where possible, figures are accurate to leading order in L. 



Circuit 


Qubits 


Depth 


Simplicity |20| 


~ 5L 


O(L^) 


Speed Ej] 


0(L 2 ) 


0(L log L) 


Qubits l22l 


~ 2L 


~ 32L 3 


Tradeoff 1 ESI 


~ 50L 


~ 2 19 L 12 


Tradeoff 2 ESI 


~ 5L 


~ 3000L 2 


LNN circuit ESI ~ 2L 


~ 32L 3 





Entanglement is a powerful resource available to quantum computers, however arbitrary entangled 
states cannot be represented efficiently on classical computers, with memory requirements scaling 
exponentially with the total number of qubits. Hence, minimising the total number of qubits is a nec- 
essary requirement for computationally tractable simulations. Also, in the short term, many current 
QC architectures face a difficult hurdle in fabricating a large number of reliable qubits, making mini- 
mal qubit circuits desirable. 

Beauregard 1221 details an implementation of the QPF subroutine appropriate for architectures al- 
lowing for the arbitrary coupling of qubits (non-LNN), in which modular addition and multiplication 
circuits are performed in Fourier space. An appropriate circuit for Linear Nearest Neighbour (LNN) 
architectures used in this investigation, detailed in Ref. |24|, uses the same method in order to re- 
duce the total number of qubits required. Both the LNN circuit and a slightly modified version of the 
Beauregard circuit require 2L + A qubits and have identical depths and gate counts to leading order in 
L. 

3 Error models and analysis 

In our simulations, errors were simulated using the discrete model in which a single qubit \<j>) = 

a\0) + (3\1) can experience a bit flip X\<f>) = a\l) + /3|0), a phase flip Z\<j>) = a\0) - /3|1), or both 
at the same time XZ\<j>) = a\l) — /3|0). These operators are simply the set {a x , ia y , a z }. 

These discrete error operators are then applied to each qubit, after each operational time step with 
probability p/3 (i.e each error has identical probability of occurrence, with the total probability of 
error given by p). The operational time for all two qubit gates is assumed to be identical and all single 
qubit gates combined with neighbouring two qubit gates via the canonical decomposition 1 25 26 27 j. 
The discrete error model represents the most common error model used within QEC analysis. This 
model oversimplifies error effects within a quantum computer in several ways. 

• The error model used is uncorrected and random. Some architectures may be more vulnerable 
to dephasing errors (Z operations), relaxation errors (X operations) or loss of qubits (this is 
particularly relevant in linear optical systems 1311 '). 

• This model does not examine the effect of systematic errors due to inaccurate gate design. 
Inaccurate two qubit gates will generally produce correlated errors over pairs of interacting 
qubits. 

• This specific error model treats memory errors and gate errors identically, which may not be 
realistic given a specific physical architecture. 
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Although this model represents a simplification of the many diverse effects that can cause errors within 
quantum computers, our interest in LNN architectures and their close adherence to this model make it 
appropriate. Furthermore, general continuous errors are equivalent to a linear combination of discrete 
errors. Correction protocols project encoded qubits onto a state that is perturbed from an error free 
state by discrete X and/or Z gates, digitising continuous errors to a discrete set. 

Using this error model, we can analytically describe the behaviour of the QPF subroutine in the pres- 
ence of severe errors. Referring to the quantum circuit used |24|, j is obtained bit-by-bit via a series 
of measurements on a master control qubit. This master qubit simulates the entire 2L qubit register 
described in section (|2ji. The QFT on this single qubit required by Eq. @]is performed through a series 
of Hadamard gates and classically controlled single qubit rotations. In a more general analysis we can 
model the entire computer as two registers, a single master qubit and the rest of the computer. 

Consider the state of the computer at a point just before the application of a controlled modular mul- 
tiplication gate. At this point the master control qubit is in an equal superposition of |0) and |1) and 
the rest of the computer is some unknown superposition, 

2 2i -l 

\4>) — — = (|0) mas ter + 1 1) master) E a m I m ) computer- (7) 
* m— 

Now apply the modular multiplication gate, which will return a new superposition state for the 
|™)computer register (when the master qubit is in the |1) state). This new superposition is denoted 
through the coefficients, {(3 m }, 

10) = -Lo) Y. a ™l m > + 7f i 1 ) E < 8 > 

V m=0 * m— 

Prior to measurement, a classically controlled rotation (9) and a second Hadamard gate is applied to 
the master control qubit. The value of 8 is dependent on the result of all previous measurements on 
this qubit. Hence the state just before measurement is, 

2 2L -1 2 2L -1 



\<t>) = \n E («™ + ei6 M\m) + i|l> Y (««• - z l6 Pm)\m). (9) 

m— m— 

With the probability of measuring a 1 or is given by, 



2 -1 

. 1 1 % 1 1 



P( 2 T 2 } = 2 ± 4 E ( ei "<A» + ^ e a m p* m ), (10) 

m— 

using, 

2 2Jj -l 2 2L -1 

Y Kn| 2 = Y \P™\ 2 = 1. (ID 
m— m— 

Errors cause the summation in Eq. ^|to asymptote to resulting in an equal probability p = (0.5) 2i 
of each j being observed. 
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The period of the function, r, dictates the number of non-zero coefficients {a m ,f3 m } and the spe- 
cific value of j simply changes the sequence of l's and O's measured at each step. Since errors act 
to randomly perturb these sets of coefficients, considering different values of r and/or j will have no 
effect on the stability of the QPF subroutine. 

The simulated QPF circuit is extremely complex and hence requires a large amount of classical simu- 
lation time. Ideally, simulations would proceed by applying a predetermined number of discrete error 
gates to every possible location within the circuit and averaging the probability of success, s, over all 
possible locations. For example, Fig. |2]shows the effect of a single X error on the QPF success prob- 
ability, s, for the first modular multiplication gate in the LNN, L = 5, circuit. From this we can see 




Fig. 2. Map showing how the location of a single bit flip error plays a major role in the final output success of 
the LNN circuit. This image is for L = 5 (14 qubits), and shows the first modular multiplication section of the 
circuit. Each horizontal block represents one of the 14 qubits while each vertical slice represents a single time 
step. Darker areas represent successively lower values for s. 

that the spacio-temporal location of an error plays a major role in the final value of s calculated, with 
various sections invariant to the bit flip error. In order to analyse the behaviour of the QPF subroutine 
we take an ensemble average over all possible error locations. For example, in fig. [2] the average 
value of s over all possible locations for a single error is s = 0.34. Most circuits are far too large to 
map out this topology efficiently: we are limited by computational resources to 50 statistical runs to 
obtain an approximate average value of s for these circuits. However, the results show that there is 
still sufficient data to observe trends in the results. 

4 Stability under a fixed number of errors 

The classical simulation algorithm employed used a state vector representation. Matrix operations 
were performed to simulate both quantum gates and error operations. In figs [5] and |4] we plot the 
success of the QPF subroutine as a function of the number of discrete errors, we plot the results for 
2L + 4 = 14, 16, 18, 20, representing factorisation of composite numbers from N = 27 to N = 247. 
Simulations examined functions that each had a period r = 6. Table[2]show the functions f(k) used 
for each value of L. These simulations aim to investigate the behaviour of the QPF subroutine for 

Table 2. Functions used for various values of L. Note that for 2L + 4 = 14, 16 the functions used are not 
products of two primes. With some slight modifications to the classical post-processing, Shor's algorithm can still 
be used to factor such numbers. Since we are only investigating the reliability of the QPF subroutine, this is not 
relevant to our analysis. 



2L + 4 


f(k) = x k modN, with r = 6 


14 


8 fc mod27 


16 


31 fe mod63 


18 


10 fc mod77 


20 


27 fc mod247 
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high component precision, close to the 1 / n p bound. Simulations were performed in a half-stochastic, 
half-deterministic manner: The type and spacio-temporal location of discrete errors occur at random, 
however we specify exactly how many errors can occur within a given run of the subroutine. 

Simulations examine the probability of obtaining the specific useful value j — [2 2L /6\. Figs [3] 
and|4]show the results for the non-LNN |22| and LNN |24| circuits respectively. For clarity, we have 
suppressed the statistical errors on these log plots. The complete data sets are given in Appendix A. 
The definition of success for the QPF subroutine, given in section |2] takes into account that many 

non — LNN QPF subroutine under fixed number of errors for various values of L. 
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Fig. 3. Plot showing the relative probability of measuring j = [2 2L /6] as a function of the specific number of 
errors for the non-LNN circuit. The curves represent L = 5 to L = 8. The horizontal lines show the point of 
random output for each successive value of L. 



different values of j may be used to determine r. However, for the sake of this analysis, we are only 
concerned if the QPF subroutine returns with high probability a value of j that is theoretically pre- 
dicted. Therefore in Figs. [3]and|4]we normalise the plots such that an error free calculation returns 
j = [2 2L /6\ with probability one and Shor's algorithm succeeds with a single call to the QPF sub- 
routine. As the number of errors increase the probability of measuring j = [2 2L /6\ decreases until it 
reaches the point of random output, at this stage the QPF subroutine performs no better than randomly 
choosing a value of j in the range j = — > j = 2 2L . 

Figs. and |4] clearly shows how the quantum speed up of the QPF subroutine, and hence Shor's 
algorithm, diminishes to a point where it is no different to randomly choosing a value from the j 
register, as the number of errors increases (represented by the horizontal lines). At this point, any 
quantum processing can no longer be identified from the probability spectrum for j . 

Accurate curve fits are extremely difficult to obtain from the limited amount of data available due 
to long computation times. Each point represents 50 separate simulations where the total number of 
errors occur randomly within the QPF circuits. In order to get sufficient data to extract meaningful 
fits for each of these curves, one would expect the number of statistical runs should be the same order 
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LNN QPF subroutine under fixed number of errors for various values of L. 
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Fig. 4. Error stability for the LNN circuit. Equivalent to fig. [3] 



as the number of possible error locations (or error combinations). For example, in the L = 5 circuit, 
for one error, the number of possible error locations and types is 18000. Hence, it is quite surpris- 
ing that even 50 statistical runs provides enough data to obtain a qualitative picture of how the QPF 
subroutine behaves for various values of L (the plots for each circuit, including statistical errors, are 
detailed in Appendix A). To reduce the statistical errors and obtain accurate curve fits for these plots, 
further simulations are required, preferably using the density matrix formalism. However, from this 
data we can still draw qualitative conclusions about the average robustness of the QPF routine as a 
function of increasing number of errors. 



To verify that a quantum computer implementing the QPF routine is processing in the quantum regime, 
it would be sufficient to observe peaks within the probability spectrum for j. The sharper the peaks, 
the fewer repetitions of QPF required and the more practical the computation. For very low visibility 
peaks, the number of repetitions of QPF scales exponentially with L, nullifying the advantages of the 
quantum algorithm over its classical version. 

These simulations show that a maximum error rate of l/n p for all problem sizes is not required to 
obtain better performance than classically searching through the j values. By inspection of figs [5] 
and an estimate can be made regarding the number of errors (as a function of L) before quantum 
processing in the QPF cannot be identified [Fig. [5). Fig. [5]represents only a preliminary estimate 
from figs. [3]and[4] additional data is required to perform an accurate curve fit. The purpose of fig. [5] 
is simply to demonstrate that when attempting to observe quantum processing, experimentally, more 
than one error can be tolerated, and the number of errors increases with L. 



When attempting to realise the full potential of the QPF routine, the probability of useful output 
should be kept as high as possible. It can be seen from figs. [3]and|4]ihat even a single error signif- 
icantly reduces this probability. Therefore, our simulations support the view that for large scale im- 
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plementation of the QPF routine, ideally no errors should occur in the circuit. This would, of course, 
be achieved though quantum error correction, with work by Steane 1 121 1291 [Tl | already examining ef- 
fective logical qubit error rates, given a specific physical error rate, for various error correcting codes. 

Maximum number of errors producing non — random output as a function of L. 
50 C - 



E 20 r 
E : 



Fig. 5. Estimate on the maximum number of errors possible for each value of L before the LNN circuit becomes 
equivalent to a classical random search. X = non-LNN circuit, [J = LNN circuit. 



The error behaviour for the LNN and non-LNN circuits are largely indistinguishable from each other. 
However, there is a slight difference in the error sensitivity of the two circuits. We attribute this to a 
minor increase in the LNN circuit depth. As expected, the overall area (n p ) of the circuit is the domi- 
nating factor in its sensitivity. The mesh circuit 1 24 1 required in the LNN design is the major difference 
between the LNN and non-LNN circuits. This section of the LNN circuit acts to slightly increase the 
overall depth, from 32L 3 + 66L 2 -2L-1 for the non-LNN circuit to 32L 3 + 80L 2 - 4L - 2 for the 
LNN design [Table. [3). Hence the sensitivity of the LNN circuit increases slightly compared with the 
non-LNN circuit. 

Table 3. Total circuit depths (K) for the LNN and non-LNN circuits, for L = 5 to L = 8 



L 


LNN Circuit 


non-LNN Circuit 


5 


5978 


5639 


6 


9766 


9275 


7 


14866 


14195 


8 


21470 


20591 



The scaling in the QPF subroutine shown by our simulations can be utilised when testing such a 
complex quantum circuit for evidence of quantum processing. As mentioned previously, peaks within 
the probability spectrum of j are indicative of quantum processing and our simulations have show that 
such peaks will be present even when component precision is not bounded by 1/ n p . 

Although we are unable fit a specific curve to the data shown in figs. [3]and|4] we can consider several 
different cases for the scaling of the tolerable number of errors. Since n p « 0(L 4 ) for the quantum 
circuit used in the simulations, and restricting our analysis to integer powers of L, we can safely bound 
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the tolerable error scaling as at most 0(L 3 ). If this scaling were 0(L 4 ) then the required component 
precision would scale as 0{L A ) / n p w constant, implying that as the quantum circuit increases size, 
the required component precision remains constant. This argument also is valid for scaling faster than 
0(L 4 ). If such a fast scaling occurred, increasing the circuit size would lead to a decrease in required 
component precision. Fig. [6]examines the required physical component precision required to observe 
quantum processing in the QPF, for potential scalings of l/n p (x), L/n p (o), L 2 jn v (A) and L 3 /n p 
(□). As a probability spectrum of j is needed in all cases, the number of total QPF routine iterations 
needed is approximately of order 2 2L . Therefore, only L = 5 to L = 10 is shown, since for higher 
L the total number of circuit iterations become prohibitively large. Table |4] examines the minimum 




Fig. 6. Required component precision in order to observe quantum processing for small values of L in the QPF 
subroutine. Each curve represents a separate type of additional scaling demonstrated by simulations, X = l/n p , 
o = L/n p , A = L' 2 /n p andQ = L 3 /n p 



physical component precision required to observe quantum processing for L = 5to_L=10 for a 
component precision of l/n p and for linear, quadratic and cubic scaling. Depending on the exact 

Table 4. Minimum component precision required to apply the QPF subroutine to at least L = 10 for various 
scalings in component precision. 



required component scaling 


component precision required 


l/np 


1.3 x icr° 


L/np 


1.3 x icr 5 


L 2 /n p 


1.3 x icr 4 


L :i /n p 


1.3 x icr 3 



nature of this scaling, quantum processing can be observed for L = 5toL=10 with a physical 
component precision between one and three orders of magnitude lower than the l/n p bound. 

Even though the computational resources required to effectively identify the functional form of this 
scaling are unavailable, we can conclude that some non-constant scaling of component precision with 
L is present leading to a reduction in the physical component precision required to observe quantum 
processing for small values of L. This result would reduce the pressure on experimental fabrication, in 
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the short term, by potentially removing the need to implement complicated error correction protocols 
to simply observe quantum processing in small instances of the QPF subroutine. 

5 Conclusion 

We have simulated the quantum part of Shor's algorithm, the quantum period finding subroutine 
(QPF). Our simulations have shown that the structure of this quantum circuit leads to a robustness 
above the naive 1 /n p approximation for component precision, if one only wishes to demonstrate the 
existence of quantum processing. Depending on the functional form of this scaling, quantum process- 
ing can still be observed for error rates significantly higher that the l/n p bound for small instances 
of the QPF subroutine. This additional robustness is advantageous in the short term since introducing 
complicated QEC protocols to simply observe quantum processing will be difficult. However, the re- 
sults of our simulations suggest that for large practical problem sizes, extensive use of error correction 
will be required to ensure error free calculations. 

The restriction to a linear nearest neighbour design does not significantly alter the sensitivity of the 
subroutine, provided that appropriate LNN circuits can be designed roughly equivalent in depth and 
qubit numbers as non-LNN circuits. 

Further work includes performing a detailed density matrix simulations in order to determine an ap- 
proximate functional form for this additional scaling, using an appropriate simulator such as QuID- 
DPro Il30l . 
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Fig. A.l. Plot showing the stability of the non-LNN (right plot) and LNN (left plot) QPF circuit for L = 5, 
included on this plot are errors associated with the stochastic nature of simulations, representing the variance on 
the 50 ran data set used in simulations. 
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n-LNN QPF subroutine under fixed number of errors for L = 6 
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Fig. A.2. Plot showing the stability of the non-LNN (right plot) and LNN (left plot) QPF circuit for L = 6, 
included on this plot are errors associated with the stochastic nature of simulations, representing the variance on 
the 50 run data set used in simulations. 
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Fig. A.3. Plot showing the stability of the non-LNN (right plot) and LNN (left plot) QPF circuit for L = 7, 
included on this plot are errors associated with the stochastic nature of simulations, representing the variance on 
the 50 ran data set used in simulations. 



LNN QPF subroutine under fixed number of errors for L=S. 
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Fig. A.4. Plot showing the stability of the non-LNN (right plot) and LNN (left plot) QPF circuit for L = 8, 
included on this plot are errors associated with the stochastic nature of simulations, representing the variance on 
the 50 run data set used in simulations. 



